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Abstract. We discuss the limiting behavior (using the notion of T-limit) 
of the 3d nonlinear elasticity for thin shells around an arbitrary smooth 2d 
surface. In particular, under the assumption that the elastic energy of defor- 
mations scales like h 4 , h being the thickness of a shell, we derive a limiting 
theory which is a generalization of the von Karman theory for plates. 
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1. Introduction 

The derivation of lower dimensional models for thin structures (such as mem- 
branes, shells, or beams) from the three-dimensional theory, has been one of the 
fundamental questions since the beginning of research in elasticity [19] . Recently, 
a novel variational approach through r '-convergence has lead to the derivation of 
a hierarchy of limiting theories. Among other features, it provides a rigorous jus- 
tification of convergence of three-dimensional minimizers to minimizers of suitable 
lower dimensional limit energies. 

In this paper we discuss shell theories arising as r-limits of higher scalings of the 
nonlinear elastic energy. Given a 2-dimensional surface S, consider a shell S h of 
mid-surface S and thickness h, and associate to its deformation u the scaled per 
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unit thickness three dimensional nonlinear elastic energy E elasUc {u,S h ). We are 
interested in the identification of the r-limit Ip of the energies: 

h- f3 E elastic {-,S h ), 

as h — > 0, for a given scaling (3 > 0. As mentioned above, this implies conver- 
gence, in a suitable sense, of minimizers u h of E elastlc (subject to applied forces) 
to minimizers of two-dimensional energy Ip, provided E elastic {u h , S h ) < Chi 3 . 

In the case when S is a subset of R 2 (i.e. a plate), such T-convergence was first 
established by LeDret and Raoult [T5] for (3 = 0, then by Friesecke, James and 
Miiller [TO] for all (3 > 2 (see also [H] for results for (3 = 2 under additional 
conditions). In the case of < (3 < 5/3, the convergence was recently obtained 
by Conti and Maggi [5], see also [2]. The regime 5/3 < (3 < 2 remains open and 
is conjectured to be relevant for the crumpling of elastic sheets. Other significant 
results for plates concern derivation of limit theories for incompressible materials 
[31 HI UH] ! for heterogeneous materials [25] and through establishing convergence of 
equilibria, rather than strict minimizers |20[ 122] . 

Much less is known in the general case when S is an arbitrary surface. The first 
result by LeDret and Raoult [16j relates to scaling (3 = and models membrane 
shells: the limit Iq depends only on the stretching and shearing produced by the 
deformation on the mid-surface S. Another study is due to Friesecke, James, Mora 
and Miiller [8], who analyzed the case (3 = 2. This scaling corresponds to a flexural 
shell model, where the only admissible deformations are those preserving the metric 
on S. The energy I2 depends then on the change of curvature produced by the 
deformation. 

All the above mentioned theories (as well as the subsequent results in this pa- 
per) should be put in contrast with a large body of literature, devoted to deriva- 
tions starting from three-dimensional linear elasticity (see Ciarlet [Tj and references 
therein). Indeed, since thin structures may undergo large rotations even under the 
action of very small forces, one cannot assume the small strain condition, on which 
the linear elasticity is based. 

The objective of this work is to discuss the limit energies for scalings (3 > 4, 
for arbitrary surfaces S. We now give a heuristic overview of our results, whose 
precise formulations will be presented in section 2. If £? e ' ostlc ( Uj s h ) ~ Ch@ , for any 
(3 > 2, one expects u to be close to a rigid motion R. This argument can be made 
precise by means of the quantitative rigidity estimate due to Friesecke, James and 
Miiller [9] (see also Lemma IQ]) . We further demonstrate that the first term in the 
expansion of u — R, in terms of h, equals an infinitesimal isometry V . That is, there 
is no first order change in the Riemannian metric of S under the displacement V. 
The corresponding bending energy, given in terms of the first order change in the 
second fundamental form of S, is the r-limit Ip if (3 > 4 (Theorem I2.3|) . This limit 
energy coincides with the so-called linearly elastic flexural shell model, derived in 
PQ from the linear elasticity theory. Our result guarantees therefore that, without 
any a priori smallness assumption on the strain, the use of the linearized flexural 
shell model is justified whenever the order of magnitude of the per unit thickness 
three-dimensional energy is hP with (3 > 4. 

When (3 = 4, also the second order in h change in the metric on S (stretching) 
contributes to the limiting energy. This change is induced by V , and additionally, by 
an "approximate second order displacement" w. This last notion involves studying 
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the finite strain space. For a similar situation where this space emerges see the 
discussion by Sanchez-Palencia [25] and Geymonat and Sanchez-Palencia 111 under 
the title of ill-inhibited surfaces, in the context of linear elasticity. In Theorems l2.ll 
and 12. 21 we derive the energy functional I4, which can be seen as a generalization of 
the von Kdrmdn theory for plates [14] , justified in terms of T-convergence in [lOj . 
Indeed, if S is a plate, then V and w are, respectively, the out-of-plane and the 
in-plane displacements (modulo a possible in-plane infinitesimal rigid motion). 

A particular class of surfaces when I4 simplifies to the bending energy is the 
hereby introduced class of approximately robust surfaces. We say that S is (ap- 
proximately) robust if any infinitesimal isometry V can be completed by a second 
order displacement to an (approximate) second order isometry. In other words, S 
can always further adjust its deformation, to compensate for the change of metric 
produced at second order. As a result, the total stretching of second order is in- 
significant and the T- limit consists only of a bending term (Theorem l2.3l) . We show 
three general examples of approximately robust surfaces: convex surfaces, surfaces 
of revolution, and developable surfaces without flat parts. An example of a not 
approximately robust surface is a plate. 

We also address the issue of external forces, depending on the reference configu- 
ration, namely the dead loads ( Theorem I2.5[) . Under a vanishing average condition 
and a suitable scaling of the forces f h applied to S , Theorem 12.11 provides an 
information of the deformation of S h assumed in response to the load. In addi- 
tion, the appropriate limit force / identifies the set of possible rotations the body 
will undergo. This phenomenon is easily observed: if f h is "compressive" , then 
S h prefers to make a large rotation rather than undergoing a compression, and an 
alignment of the infinitesimal isometry V with the force is energetically preferable. 

The identification of T-limit for any scaling in the range j3 £ (2, 4) and non-flat S 
is still open. In analogy with the analysis developed in [TO] for plates, the construc- 
tion of a recovery sequence requires finding an exact isometry of S, coinciding with 
a given second order isometry. Another direction of study concerns shells, whose 
mid-surface is inhibited (or infinitesimally rigid). Examples of such are closed or 
partially clamped elliptic surfaces. In this case the limit functionals that our theory 
yields are identically equal to zero. This suggests looking for higher order terms 
in the development of the three-dimensional energy in the sense of T-convergence. 
These are subtle issues and we plan to address them in a forthcoming paper. 

Acknowledgments. We thank Stefan Muller for helpful discussions. A large 
part of this work was carried out while the second author was visiting the Insti- 
tute for Mathematics and its Applications in Minneapolis (USA), whose support 
is gratefully acknowledged. M.L. was partially supported by the NSF grant DMS- 
0707275. M.G.M. was partially supported by the Italian Ministry of University 
and Research through the project "Variational problems with multiple scales" 2006 
and by GNAMPA through the project "Problem! di riduzione di dimensione per 
strutture elastiche sottili" 2008. 

2. An overview of the main results 

Let S be a 2-dimensional surface embedded in M 3 . We assume that S is com- 
pact, connected, oriented, and of class C ,x , and that its boundary dS is the union 
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of finitely many (possibly none) Lipschitz continuous curves. Consider a family 
{S h }h>o of thin shells of thickness ft around S: 

S h = {z = x + tn(x); x e S, -ft/2 < t < ft/2}. 

We will use the following notation: n(x) for the unit normal, T X S for the tangent 
space, and U(x) = Vn(x) for the shape operator on S, at a given x £ S. The 
projection onto S along n will be denoted by n, so that: 

tt(z)=x Vz = x + tn(x) e S"\ 

We will assume that ft < ft , with ft > sufBciently small to have 7r well defined, 
and so that: 1/2 < |Id + ffl(x)| < 3/2 for all z as above. 

For a W 1,2 deformation of a thin shell u h : S h — > M 3 , we assume that its elastic 
energy (scaled per unit thickness) is given by the nonlinear functional: 

E elastlc (u h ,S h ) = \ ( W{Vu h ), 
ft Js h 

where the stored-energy density function W is nonnegative and C 2 in some open 
neighborhood O of SO(3), in the space M 3x3 of 3 x 3 real matrices. Moreover, W 
is assumed to satisfy, for all F € R 3x3 and some C > 0: 

W(RF)=W(F) Vi?eSO(3), 

W(R)=0 VReSO(3), 

W{F) > Cdist 2 (F,SO(3)). 

Here SO(3) denotes the group of proper rotations. Recall that the tangent space 
to 5*0(3) at Id is the space of skew-symmetric matrices: 

so(3) = {F e M 3x3 ; F = -F T }. 

It is convenient to view u h through their rescalings y h E W /1 ' 2 (S I?! ' ,R 3 ): 

y h (x + tn(x)) = u h (x + th/h a n(x)) Vx e S Vt e {-h /2, h Q /2). 

The advantage is that all y h have the common domain S h ° . We are concerned with 
the limiting behavior of the energies: 

i h (y h ) = \l w(Wu h ), 

h Js h 

relative to low energy deformations. That is, we want to discuss the limit, as 
ft — > 0, of the functional I h /e h 7 for a given sequence of positive numbers e h , which 
we assume to satisfy: 

(2.1) lim e' l /ft 4 = k 2 < oo. 

h — ^0 

With this in mind, define the related scaled average displacement: 
(V h [y h ])(x) = -JL= -f ' 2 y h (x + tn)-x dt. 

Ve" J -h a /2 

Since we will frequently deal with such vector fields V € W 1 ' 2 (S, M 3 ) on the surface, 
we introduce the following notation. By sym VV(x) we mean a bilinear form on 
T X S given by: (sym VV{x)t)j] = ±[(d T V(x))ri + [d v V{x)) T ], for all t,t] e T X S. 
Also, given a matrix field A e L 2 (S, M 3x3 ), by A tan (x) we denote the tangential 
minor of A at x € S, that is [(A(x)T)r]] TiTie T x s- 
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Our first main result is the following: 

Theorem 2.1. Assume &2.1\) and let u h € W 1,2 {S h , R 3 ) be a sequence of defor- 
mations such that the sequence of scaled energies 

{-prl h {y h )} is bounded. Then 
there exist rigid motions of Mr, given through proper rotations Q h e 50(3) ana 
translations c h G R 3 such that for the normalized deformations: 

y h {x + tn) = (Q h fy h (x + tn) - c h 

the following holds. 

(i) y h converge in W 1,2 (S ha ) to tt. 

(ii) V h [y h ] converge (up to a subsequence) in W 1 ' 2 (S) to some vector field V G 
W 2,2 (S, R 3 ) with skew- symmetric gradient, that is: 

(2.2) d T V(x) = A(x)t Vt g T x S, a.e. x£ S 
for some matrix field A G W 1 ' 2 (S, R 3x3 ) such that: 

(2.3) A(x) G so(3) Va; G S. 

(iii) j^sym W h [y h ] converge (up to a subsequence) weakly in L 2 (S) to some 
symmetric matrix field B ta n on S. 

(iv) There holds: 




where: 

(2.4) 

I{V, B tan ) = \j s Qt (x, B tan - |(A 2 ) tan ) +1^J^Q 2 (x, (X7(An) - AW) tan ) . 

The following quadratic, nondcgcncratc forms are of relevance here: 

(2.5) Q 3 (F) = D 2 W(ld)(F,F), Q 2 {x,F tan ) = min{Q 3 (F); (F — F) tan = 0}. 

The form Q3 is defined for F G R 3x3 , while Q.2(x,-), for a given x G S is defined 
on tangential minors F tan of such matrices. Both forms Q3 and all Q,2{x, •) are 
positive definite and depend only on the symmetric parts of their arguments (see 
0). 

Theorem 12.11 will be proved in sections [3] and [4] One of the crucial ingredients 
is a result on approximating large deformations [S]. For completeness, we sketch 
its proof, in the setting of shells, in Appendix A. We also note that because of 
the non trivial geometry of the shell, the limiting energy /, in general, exhibits a 
dependence on the point x G S, although the three-dimensional configuration is 
homogeneous. 

Our second main result concerns the possibility of recovering the functional 
I(V,B ta n) in (|2.4p (or its components), as the limit of scaled energies -prl h (y h ), 
for some sequence of deformations. For this, define the finite strain displacement 
space: 

B = { limsym Vw h ; w h G W 1,2 (S, R 3 )|, 

where limits are taken in L 2 (S) (clearly, both the weak and the strong convergences 
yield the same 23). We then have: 
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Theorem 2.2. Assume 112.1)) . For every V € W 2 > 2 (S, R 3 ) satisfying i2.2]) and 
12. 3\) . and every Btan € B, there exists a sequence of deformations u h £ W 1,2 (S h , R 3 ) 
such that: 

(i) y h converge in W x ' 2 {S h °) to it. 

(ii) V h [y h ] converge in W^ 2 (S) to V. 

(iii) Tjsym W h [y h ] converge in L 2 (S) to B tan . 

(iv) Recalling the definition \2.J$ one has: 

limJ K I h (y h ) = I(V,B tan ). 
h^o e n 

The form of the limiting energy functional / simplifies, when the space B is large 
enough to choose B tan so that the first term in (|2.4|) vanishes. That is, we call S 
"approximately robust" if for every V £ W 2 ' 2 (S,R 3 ) satisfying ([2~2]l (f2~3]) . one has 
(A 2 )tan S B. Then we have: 

Theorem 2.3. Assume 12.1]) . Let k = or let S be approximately robust. Then 
for every V € W 2 ' 2 (S, R 3 ) satisfying &2.2\) and &2.3\) , there exists a sequence of 
deformations u h £ W 1 ' 2 (S h , R 3 ) such that (i) and (ii) of Theorem [Ql hold. More- 
over: 

limi/V) = /(n 

where 

( 2 -6) I(V) = ^J s Q 2 (x,iy(An)-AU) tan ). 

Theorems 12.21 and 12.31 will be proved in section [5J In section O we discuss the 
space B and approximately robust surfaces. In particular, we shall see that convex 
surfaces, surfaces of revolution, and non-flat developable surfaces are approximately 
robust. 

Theorems 12.11 and 12.21 (or 12. 3[) can be summarized (although they provide more 
information than the below statement), using the language of T-convergence. For 
completeness, the following result will be presented in Appendix B. 

Corollary 2.4. Assume \2.1\) . 

(i) Define a sequence of functionals: 

T h : W h2 (S h °,R 3 ) x W 1 ' 2 (S,R 3 ) x L 2 (S, R 2x2 ) — ►! 

T h (y h , V\Bl n ) = I -jJ h (v K ) V Vh = and <n = bv m W ' l > 

\ +oo otherwise. 

Then T h Y -converge, as h — > 0, to the following functional: 

( I(V, B tan ) if V = ?r, V £ W 2 - 2 satisfies [KM), fOU . 
T(y, V, Btan) = I and B tan £ B, 

I +oo otherwise. 

(ii) Assume that k — or let S be approximately robust. Define the functionals: 

T h : W 1 - 2 (S ho ,m 3 ) x W 1 ' 2 {S,R 3 ) — >R 

T\y\V h ) = \ h lh{yh) lf Vh = Vh ^ 
I +oo otherwise. 
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Then T h Y -converge, as h — > 0, to the functional: 

y\ = { I(V) if y = TT andV E W 2 - 2 satisfies (EM, WM, 
\ +00 otherwise. 

All statements above remain valid if the product spaces (the domains of functionals 
T , T ) are equipped with the weak (instead of strong) topology. 

We further consider a sequence of forces f h S L 2 (S h ,lS?), acting on thin shells 
S . For simplicity, we assume that: 

f h (x + tn(x)) = hVeF det (Id + <n(x)) _1 /(x), 

where / € L 2 (S, R 3 ) is normalized so that: 

(2.7) f f = 0. 

Js 

Define m to be the maximized action of force / on S over all rotations of S, and 
let M. be the corresponding set of maximizers: 

(2.8) M = \q e SO(3); [ f(x) ■ Qx dx = m = max [ f ■ Qx 

I Js QeSO(3)J s 

The total energy functional on S h is given through: 

J h (V h )=l\y h )+m h -~ j f h u\ 



S h 

where m h = hVe^m. 

Theorem 2.5. Assume K2l\) and (E7\ ). Then: 

(i) For every small h > one has: 

> inf j-^V 1 ); u h S W^ 2 (S h , R 3 )J > -C. 

(ii) If u h E W 1 ' 2 (S, R 3 ) is a minimizing sequence of -jrJ h , that is: 

(2.9) 

i/ien i/iere exists Q h E SO(3) and c h E R 3 swc/i i/ia£ /or i/ie normalized 
deformations y h — (Q h ) T y h — c h the convergences of Theorem \2.1\ (i) (ii) 
and (Hi) hold. The convergence of (a subsequence of) j^symS7V h [y h ] to 
Btan in (Hi) is strong in L 2 (S). 

Moreover, the set of accumulation points of {Q h } is contained within A4. 
Any limit (V, B ta n,Q) minimizes the functional: 



J(V, B tan , Q) = I(V, B tan ) - [ f-QV, 

Js 



over all V E W 2 ' 2 (S,R 3 ) satisfying (EM (EM, all B tan eB andQ E M. 
(iii) If k — in (EM, 0T if S is approximately robust, then for any minimiz- 
ing sequence as in (EM, we obtain convergences of y h , V h [y h ] and Q h as 
described in (ii) above, and the limit (V,Q) minimizes the functional: 



J(V,Q)=I(V)- f-QV 
Js 

over all V E W 2,2 (S, M 3 ) satisfying (EM WM> and all Q E M. 



s 
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In section we prove Theorem 12.51 and explain the significance of the set M. in 
the setting of dead loads. 

The lower bound on the functionals J and J, as well as attainment of their 
infima, can be proved independently, under conditions (|2.7p and: 



(2.10) / f(x) ■ QFx dx = VQeM VF G so(3). 

Js 

Here M. is any closed, nonempty subset of 50(3). When M. has the form (|2.8p . 
then (|2.10p follows from (|2.8|1 and can be seen as its linearization. This analysis 
will be carried out in Appendix C. 

3. Convergence of low energy deformations 

In this section we derive some bounds on families of vector mappings {u h }h>o, 
defined on S h , under the assumption of smallness on their energy. In what follows, 
by C we denote an arbitrary positive constant, depending on the geometry of S but 
not on h or the vector mapping under consideration. In all proofs, the convergences 
are understood up to a subsequence, unless stated otherwise. 

We will work under the following hypothesis: 



(H) 



A sequence of vector mappings u h g W 1,2 {S h , K 3 ) and a sequence of 
positive numbers e h satisfy, for small h > 0: 

(i) ~ [ W(Vu h ) < Ce h , 
h Js" 

(ii) lim e h /h 2 = 0. 



As for the flat case in [TO], the first crucial step is the following approximation 
result: 

Lemma 3.1. For each u h as in (H) there exist a matrix field R h € W 1,2 (S, K 3x3 ) 
such that: 

R h (x) e so(3) Vx e s, 

and a matrix Q h £ SO(3) such that: 

(i) \\Wu h - R h ir\\ L 2 {sh) < Ch 1 ' 2 ^, 

(ii) \\VR h \\ms) ^Ch' 1 ^, _ 

(iii) ||(g' i ) T i?' l -Id|j LP( s) <Ch- x ^, for allpe [l,oo), 

The proof follows from Lemma |8. II given in Appendix A, in view of: 

E{u'\S h )= ( dist 2 {X7u h ,SO{3))<C [ W{\7u h )<Che h 
Js h Js h 

so that lim/^o h^E(u h ,S h ) = by hypothesis (H). 

Lemma 3.2. Assume (H) and let R h ,Q h be given as in Lemma \3.1\ There holds: 

(i) lim (Q h ) T R h = Id, m W 1 - 2 ^) and in L*>(S). 

Moreover, there exists a W 1,2 skew- symmetric matrix fields A : S — > so(3) such 
that: 

(ii) lim ({Q h ) T R h - Id) = A, weakly in W 1 ' 2 (S) and (strongly) inD>{S). 
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(iii) lim ^sym i(Q h ) T R h - Id) = \a 2 , in D>(S). 
h^o e n 2 

In (ii) and (iii) convergence is up to a subsequence (that we do not relabel). In (i), 

(ii), and (iii) the appropriate convergence holds for all p £ [1, oo). 

Proof. The convergences in (i) follow from Lemma 13.11 in view of (H). To prove 
(ii), notice that the sequence: 

A h = -JL ((Q h ) T R h - Id) 



is bounded in W 1,2 (S) and so it has a weakly converging subsequence. By compact 
embedding of W 1,2 (S) into L P (S) the convergence is strong in L P (S). One has: 

A h + (A h f = -A= ((Q h ) T R h + (R h ) T Q h - 2Id) = -^-(A h f ■ A h . 



The latter converges to in L P (S), and therefore the limit matrix field A is skew- 
symmetric. The above equality proves as well that: 

lim — |L= sym A h = -A 2 

in L P (S), which implies (iii). H 
Recall the rescaling: 

y h (x + tn(x)) = u h (x + th/h n(x)) Va; S S Vi G (-h /2, h /2), 
so that y h £ W 1 > 2 (S h °,M 3 ). Also, define: 

V h y h (x + tn(x)) = Vu h (x + th/h n(x)) 
By a straightforward calculation we obtain: 
Proposition 3.3. For each x S S, t S (— ho/2, Iiq/2) and r £ T X S there hold: 
d T y h (x + tn) = V h y h (x + tn) (Id + th/h U(x)) (Id + tll(x))" 



dny h {x + tn) = j-Vhy h (x + tn) n(x). 

Moreover, for I h {y h ) = ^ J sh W(\7u h ) one has: 

I h {y h ) = t I W{X7 h y h (x + tn)) • det [(Id + th/h U) (Id + tn)" 1 ] 
ho Js h o 

= f W{V h y h (x + tn)) ■ dct [Id + th/h U(x)} dt dx. 

Js J-h /2 

Also, directly from Lemma \3. II (i) and Lemma [3T2l fii) there follows: 
Proposition 3.4. Assume (H). Then: 

(i) \\V h y h - R\\\ LHS u 0) < CVe*. 
h 

' e 



(ii) lim —= {{Q h ) T ^hV h ~ Id) = An, in L 2 (S h °) up to a subsequence. 
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We will consider the corrected by rigid motions deformations y h <E W 1 ' 2 (S h °, 
and averag ed displacements V h £ W^&R 3 ): 



y h = (Q h fy h -c h 



V 



V h [y h ] = 



h /2 



y h (x + tn) - x dt, 



ho/2 



where c^f s /_^ 2 /2 (Q» 



dt dx, so that f s V h = 0. 



Lemma 3.5. As 



(H). Then: 



(i) lim y h = tt, in W 1 ' 2 (S h °). 

h—>0 

(ii) lim V h = V, in W^ 2 {S) up to a subsequence. 

h—*Q 

The vector field V in (ii) has regularity W 2 ' 2 (S, R 3 ) and it satisfies d T V(x) — A(x)t 
for all t £ T X S . The W 1,2 skew- symmetric matrix field A : S — > so(3) is as in 
Lemma \3.S\ 



Proof. 1. In view of Proposition 13.31 and Proposition 13.41 we have 
1 Vteny" - ((Q h ) T R h ) tan ■ (Id + th/h U)(ld + tU) 
J^ll^o) < Ch\\V h y h \\ LHSh0) < Ch. 
To prove convergence of V h , consider: 



(3.1) 



\\L 2 (S h o) - ° v c 



VV h {x) 



h r ho/2 

-F=f V tan y h {x + tn)(\d + m) -Iddt 
ve" J -ho/2 



(3.2) 



'-ho/ 
rho/2 



■f ° Ntanf ~ ((Q k ) T R h ) tan (M + ffl)- 1 ) (Id + til) df 
J-ho/2 K ' 



-ho/2 

{Q h ) T R h {x)-ld 



We see that by (13. 1 [) the first term in the right hand side above converges to in 
L 2 {S h °), as h — * 0. The second term converges, up to a subsequence, to A tan by 
Lemma T3. 21 (ii). Therefore W h converges to A tan in L 2 (S) and since j s V h = 0, 
we may use Poincare inequality on S to deduce (ii). 

2. To prove (i), notice that by (|3.ip and Lemma [3.21 we obtain the following 
convergences in 1? (S h ° ) : 

lim V tan y h = (Id + ill)- 1 = \7 tan TT, 
lim dnf = 0. 

h^0 

Therefore \7y h converges to Vir in L 2 (S h °). 

Since the sequence {V h } is bounded in L 2 (S), it also follows that: 

rho/2 



(3.3) 



lim 

7i->0 



ho/2 



y h - 7T dt 



= 0. 



L 2 (S) 



Now, let g(x + tn) = |det (Id + ffl(x))| 1 . Consider the two terms in the right hand 
side of: 



y - 7r lli=(s"o) 



< 



.9 



(y - tt) 



L 2 (S h o) 



.9 



(r - tt) 

/ s h a 



NONLINEAR. SHELL THEORIES 



11 



The first term can be bounded by means of the weighted Poincare inequality, by 
||V(y' 1 — 7r) 1 1 £2(5^0) and therefore it converges to as h — ► 0. The second term 
converges to as well, in view of Q3.3|) and: 







r pho/2 








/ / y h - dt dx 


< c 


JS h o 




Js J-ho/2 





h /2 



y h - 7T dt 



h /2 



L 2 (S) 



This justifies convergence of y to tt in L (S °) and ends the proof of (i). ■ 

Towards the proof of Theorem 12.11 we need to consider the following sequence 
of matrix fields on S h ° : 



G = 



1 



(R h ) T V h y n - Id 



In view of Proposition 13.41 (i), 2symG' 1 is the Ve^ order term in the expansion of 
the nonlinear strain (Vit ) T Vw , at Id. This expression will also play a major role 
in the expansion of the energy density at Id: W(Vhy h ) — W(ld + Ve h ~G h ). 

Lemma 3.6. Assume (H). Then the sequence {G h } as above has a subsequence, 
converging weakly in L 2 (S °) to a matrix field G. The tangential minor of G is, 
moreover, affine in the n direction. More precisely: 



Vr € T X S 



G(x + tn)r = G (x)t + — ■ h(An)(x) - AU{x)) r, 
ho V / 



where Gq(x) 



ho/2 



G{x + tn) dt. 



-ho/2 



Proof. 1. The sequence {G h } is bounded in L 2 (S h °) by Proposition l3.4l (i). There- 
fore it has a subsequence (which we do not relabel) converging weakly to some 
G. 

For a fixed s > 0, consider now the sequence of vector fields f s ' h G W 1,2 (S h ° ,R 3 ): 
f s ' h (x+tn) = —j= (/ioy^(^+(i+s)n)-/i(x+(t+.s)n)^-^/i t/' l (x+tr?)-/i(a;+in)^ 



We claim that f s,h converges in L 2 (S h °) (up to a subsequence) to (Ari)ir as h — > 0. 
Indeed, using Proposition 13.31 one has: 



f s ' h (x + tn) = -L= j- (h dfiy h (x + an) - hnj da 

h = £ +S {{Q h ) T V h y h {x + art) - Id) ft da, 



and the convergence follows by Proposition 13.41 (ii) . 

2. We claim that this convergence is actually weak in W 1,2 (S h °). First, notice 
that the normal derivatives converge to in L 2 (S h °) by Proposition 13.41 (ii): 



On 



H f s ' h (x + tn) = -A= (Q h ) T h h y h (x + (t + s)n) - V h y h {x + tn))n(x). 
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We now find the weak limit of the tangential gradients of f 3 ' . By Proposition ^. 31 
there holds, for all r £ T X S: 

d T f sJ \x + tfi) = — ^= (h X7y h (x + (t + s)n)(Id + (t + s)n)(Id + til)- 1 
sVe h v 

- h Wy h (x + tn) - hsUild + my^T 
= -^= (Q h f hhV h {x + (t + s)n) - V h y h (x + tn)) (Id + th/hott)(Id + m)~ l r 
+ —= ((Q h ) T V h y h (x + (t + s)n) - Id)sII(Id + tH)~ 1 T. 
By Proposition [313] (ii), the second term in the right hand side above: 

-j= ({Q h ) T v h y h (x + (t + s )n) - id)n(id + tn)- 1 

converges in L 2 (S h °) to AII(Id + tll)~ 1 . 

On the other hand, the first term equals to: 



s 



(Q h ) T R h (G h {x +(t + s)n) - G h (x + in)) (Id + th/h tt)(Id + til)- 1 
and by Lemma 13.21 (i) it converges weakly in L 2 (S ha ) to 

^ (g(x + (t + s)n) - G(x + tn)) (Id + tn)" 1 . 



ho 

s 

This establishes the (weak) convergence of f s ' h in W 1 ' 2 (S h °) 



3. Equating the weak limits of tangential derivatives, we obtain, for every r £ 
T X S: 



d T {An){x) ■ (Id + nTT 1 = — (G(x + (t + s)n) - G(x + tn)) (Id + ffl^T 

+ ^n(id + tn)- 1 r. 

This proves the lemma. I 
Finally, we have the following bound for convergence of the scaled energies J . 
Lemma 3.7. Assume (H). Then: 

]immf±-I h {y h )>~ f Q 2 (z, (sym G ) to „) + ^ / Q 2 (x, (V(Aft) - AU) tan ) . 
h^O e n 2 J s 24 J s 

Proof. By the frame invariance property of W we have: 

W{V h y h ) = W{{R h ) T V h y h ) = W(ld + Ve T G h ). 



Consider the sets D, h = {x £ S h °; (e h )^ 4 \G h (x)\ < 1}. Clearly the sequence of 
characteristic functions \n h converges to 1 in L 1 (S ,/l °), as {(e' i ) 1 / 4 G' i } converges 
pointwise to 0. Since W is C 2 in a neighborhood of Id, then by the above calculation, 
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in ilh (for h sufficiently small) there holds: 

\w{V h y h ) = ~D 2 W(Id)(Ve*G h ,Ve*G h ) 



e h 2 

el 



(3.4) 



(1 - s) \D 2 W(Id + sVe^G h ) - D 2 W{ld)] ds (G h ,G h ) 



lQ 3 (G h )+o(l)\G h \ 2 . 



Above o(l) is the Landau symbol denoting any quantity uniformly converging to 
0, as h — > 0. In view of Proposition 13 . 31 we now obtain: 

rh /2 



liminf -rl h (y h ) > liminf 1 



, , , xn h W(V h y h )det[ld + th/h U]dtdx 
e" J s J_ ho/2 

° Xn h \w{V h y h )dtdx 

S J -ho/2 e 
h /2 



lim inf / -r 
JsJ- 

liminf-/-/- Q 3 (sym( X n h G h ))+o(l) / |G h | 

h->0 2 J s J -ho/2 J S h 



1 

> - 

- 2 



ho/2 



Q 3 (sym G). 

/S J -ho/2 

The last inequality follows from positive definiteness of Q3 on symmetric matrices, 
and the fact that xo. h G h converges weakly to G, in L 2 (S h ° ) . 
By the definition of Q2 and by Lemma 13.61 we get: 

rho/2 1 r fho/2 



7: / f Q 3 (symG) = - / f Q 2 ( I ,(symG) tan ) 

Z JS J -h /2 Z J -ho/2 

r fho/2 fho/2 .2 

f Q 2 (a;, (sym G ) to „) + / f - 2 -Q 2 {x,{V{An) - AIl) t 

Js J -h /2 JS J -h /2 % 



which proves the result. H 

4. A proof of Theorem 12.11 and some explanations 

To complete the proof of Theorem 12. li in view of Lemma 13.51 and Lemma 13.71 
it remains to understand the structure of the admissible matrices Go- This is the 
content of the next lemma. 

In addition to the hypothesis (H), we now also assume the existence of the finite 
limit: 



(4.1) 



lim Ve^/h 2 < oo. 



h^O 



When e h w hP ', this corresponds to the case (3 > 4, with k > for (3 = 4 and n = 
for (3 > 4. 

Lemma 4.1. Assume (H) and Lei Go be the matrix field on S, as in Lemma 

Then we have the following convergence, up to a subsequence, weakly in L 2 (S): 

lim isym W h = (sym G + -A 2 ) 

h^Oh V 2 J tan 

where the subscript tan denotes, as usual, the tangential minor of a given matrix 
field on S. 



(4.2) 
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Proof. We use the formula (13. 2p to calculate ^sym W h . The last term in the right 
hand side gives: 

1 sym {{Q h ) T R h Id) t<w = W ^ " Id Ln - 



which converges in L 2 (S) to «;/2(A 2 ) ta „ by Lemma EJ (iii) . 

To treat the first term in the right hand side of (|3 . 2(1 . notice that for every 
t e T X S: 



1 



h /2 



h /2 



Vy n (x + tn){ld + m) - {Q n ) I R fl (x) dt 



h\T 



ho/2 



V h y h (x + tn)- R h {x) dt + 



-h /2 



h /2 



th/h Q W h y h U dt 



-ho/2 



(Q h ) T R h (x) 



fho/2 

\ (R h ) T V h y h - Id dt 

J-ho/2 



h/h 



h\T 



(Q h ) 



fho/2 

f t(V h y h -R h n) dt 

J-ho/2 



IL(x)r, 



where we used Proposition [231 Now, the second term in the right hand side above 
converges in L 2 (S) to 0, by Proposition 13.41 (i). Further, the matrix in the first 
term equals to: 



(Q^R^x) 



ho/2 



G h {x + tn) dt, 



and by Lemma 13.21 (i) and Lemma 13.61 it converges weakly in L 2 (S) to G . This 
completes the proof. ■ 

We now comment on the regularity and role of various quantities containing V 
and A, intrinsically related to the geometry of the problem. 

Remark 4.2. Notice first that if a vector field V G W 2 ' 2 (S, K 3 ) has skew-symmetric 
gradient: 

r • d T V{x) = Vt G T x S, 
then it uniquely determines a W 1 ' 2 matrix field A : 5* — ► so(3) by: 
Vr G T x S At = d T V, 

An = TL-V tan -Vt<m(Vn). 



(4.3) 



Regarding the regularity, write V as the sum of its tangential and normal com- 
ponents, to obtain: 

V = Vtan + {Vn)n, sym VV = sym VV tan + (Vn)U. 

Hence, assuming sufficient regularity of S (say, S is C 3,1 up to its boundary) it 
follows that sym W tan = — (Vri)II G W 2 - 2 (S, R 2x2 ). Using the same calculations 
as in [T] page 119, we may deduce that the tangential component V tan enjoys higher 
regularity than the vector field V. Namely, V tan G W 3 ' 2 (S, K 3 ) and: 

||Vton||w 3 . 2 (S) < C^ll^onllw 1 ' 2 ^ + II ^11 IV 2 . 2 (S) 
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By Korn's inequality, one can replace the W 12 norm of Vtan by a term of the order 
HVtanlU 2 + ||sym VVtan||i 2 , so that we finally obtain: 

1 1 Vtan 1 1 W 3 > 2 (S) ^ C(||^on|U 2 (5) + II Vft\\ W 2 - 2 (S)J ■ 

For an elementary derivation of Korn's inequality on S from Korn's inequality on 
open sets, see e.g. [17] . 

In the same manner, one can prove the following useful bound, valid under C 2 ' 1 
regularity of S : 

(4-4) \\V t an\\w 2 - 2 {S) < c[\\V t an\\L 2 (S) + 1 1 Vft\ \ w i,2 (s) \ . 

Remark 4.3. 1. The (scaled) t - derivative of Gr, which is also the argument of 
the second term in the definition of / (and /) , may be written as: 

(V(An) - AU)t = [(V(An) - AU)t] = (d T A)n. 

J tan 

This expression measures the difference of order h between the shape operator II 
on S and the shape operator II of the deformed surface Sh — (id + hV)(S) (see 
Figure rj~Tj) . To see this, let x £ S and let ti, t 2 £ T X S be such that n(x) = n X t 2 . 
The tangent map of the deformation (j> h (x) — x + hV(x) equals Id + hA, and we 
obtain the following expansion of the (scaled) normal vector n h to Sh at the point 
(f> h (x): 

n h = (d T1 4> h x d T2 4> h ^j (x) = n{x) + h{r x x d T2 V + d T1 V x r 2 ) + 0{h 2 ) 
= n + hAn + 0(h 2 ), 

where we used the Jacobi identity for vector product and the fact that A £ so(3). 
Note that \n h \ = 1 + 0{h 2 ) and therefore 

II h (Id + hA)r = d T Cr^J) = d rn h + 0(h 2 ). 

Hence the amount of bending of S, in the direction of r £ T X S, can be estimated 
by: 

(Id + /i.^)- 1 n' l (Id + hA)r -TJt = (Id + hA)~ l (d T n h + 0{h 2 )) - Ut 



= (Id + hA)-' ( (Id + hA)nr + h{d T A)n + 0{h z )j - Ht 
= (Id - hA)h(d T A)n + 0(h 2 ) = h{d T A)n + 0{h 2 ). 
A closely related heuristics is the following. By Proposition [33] (for simplicity, we 



assume here that e = h ) the tangent map Vu (i), at x £ S, is approximately a 
rotation R h (x) £ 50(3). Hence, ri h ss R h n. Assuming that limQ' 1 = Id, we may 
think that R h (x) ss Id + hA(x). The difference of the shape operators on Sh and S 
satisfies: 

{R h ) T Vn h -n m (Id + hA T ){li + hV tan {An)^j - n 

« hVtan(An) + hA T U = h(v(An) - AU 



2. In turn, the role of the first term in the definition of /: 

1 K 
(sym G ) tan = lim -sym VV h ~ - (A 2 ) ta 
h— >o n Z 
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Figure 4.1. The mid-surface S and its deformation. 



is to measure the difference of order h 2 between the metric on S and the metric 
of the deformed mid-surface. Notice that under the deformation id + hV , as in 
Figure 14. 1[ there is no first order change in the length of curves on S because the 
gradient held \7V is skew-symmetric. In geometrical terms, vector fields V with 
this property are known as infinitesimal isometries (see [27] . Chapter 12). 

Under the same condition (for simplicity we again assume that e h — ft 4 ), the 
amount of stretching of S, in the direction r £ T X S and induced by the deformation 
4> h = id + hV + h 2 w has indeed the following expansion: 

\d T (f> h \ 2 - \t\ 2 = h 2 (2tO t w + \d T V\ 2 ) + 0(h 3 ) 

= 2/i 2 (V T (symVu;)T - \t t A 2 t) + 0(h 3 ). 

5. The space of finite strains B and three examples of approximately 

ROBUST SURFACES 

The space of limits as in the left hand side of (|4.2jl plays an important role in 
defining the exact limiting energy functional on S. With this in mind, we introduce: 

Definition 5.1. The space of finite strains is the following closed subspace of L 2 (S): 
B = | lim sym Vw h ; w h € W 1,2 (S,R 3 )| 

where limits are taken in L 2 (S). 

Clearly, by Mazur's theorem, B contains all weak L 2 (S) limits of symmetric 
gradients of W 1,2 vector fields on S. 

As we shall see in Theorem I2.3[ the form of the limiting energy functional sim- 
plifies, for surfaces with large space B. 

Definition 5.2. We say that S is approximately robust, if for every V € W 2,2 (S, R 3 ) 
satisfying H2.S\) and 112. 3\) , one has: (A 2 ) tan € B. 

According to our terminology, S would be called "robust" if every admissible 
(A 2 )tan as above, equaled sym Vw for some w S W 1 ' 2 (S, M 3 ). The notion of robust 
surfaces will arise at lower scalings, that is when n = oo, which we do not consider 
in this paper. 
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Remark 5.3. An equivalent construction of B is the following. Define the linear 
space of finite strain displacements: 

W = W h2 {S,R 3 )/{w e W 1 ' 2 ; sym Vw = 0}. 

It can be normed by ||[u>]||w = ll s y m V w IIl 2 (s)- Then (B,\\ ■ \\l 2 {s)) is linearly 
isometric to the completion W of (W, || • \\w), so that the elements of B can be seen 
as generalized symmetric gradients of elements of W. 

Such construction is used in 11] or in 1 in the context of derivation of membrane 
theories from linear elasticity. Note the different regularity of the kernels considered 
in [TT] and the related explanations in [TJ, page 262. 

Remark 5.4. In general, it is complicated to directly determine the exact form of 
8 or W. The crucial step in identifying W is finding the optimal norm || • || for 
which a Korn-Poincare type inequality 

(5.1) inf U\u - w\\ ; w e W 1 ' 2 , sym Vw = oj < C||symVu||i2(s) 

holds with a uniform constant C, for all u 6 W 1 ' 2 (S, M 3 ). Unlike in the case of 
tangent vector fields, this optimal norm is usually weaker than L 2 . The reason is 
that the boundedness of the left hand side in: 

sym Vw h = sym Vu>f ajl + (w h n)Il 

does not, in general, imply L 2 boundedness of both terms in the right hand side. 

This is the case, for example, when S is (a piece of) a cylinder S 1 x [—1/2, 1/2]. 
Let ri and t 2 be the tangent unit vector fields, respectively orthogonal and parallel 
to the axis X3 of the cylinder. One can show that there exists a sequence [w h ] S W 
converging in W, such that for any choice of representatives w h £ W 1 ' 2 (S,M 3 ), 
the norms 7"i ||z, 2 (S') an d || w h n\\ h~ 1 (S) blow up. However, this is the worst case 
scenario, and one has: 

W= {v e £>'(S*,R 3 ); vti E H-^S), vr 2 eL 2 {S) 7 vneH- 2 {S) 7 
sym Vv e L 2 (S), 

1/2 pl/2 

X3(vri) = const, / vn = const, 
-1/2 i-1/2 



vt 2 = / x 3 (vn) = oj. 



VTl ~ 

In this particular case, however, as we will see below, W is isometric to the space 
of all symmetric L 2 matrix fields Btan on S. 

Remark 5.5. Flat surfaces Set 2 are not approximately robust. Indeed: 

B = {sym Vw; w € W 1 ' 2 (S, M 3 )} 

= {Btan £ L 2 {S,R 2x2 ); B T tan = B tan , curFcurl B tan - o}. 

On the other hand, given V € W 2 < 2 (S,R 3 ) satisfying (|2T2")) and (TQ|) . one has 
(A 2 )tan € B if and only if V 3 = Vn solves the degenerate Monge- Ampere equation: 
detV 2 (V 3 ) = (see [TO]). 
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A particular class of surfaces which are approximately robust are surfaces for 
which: 

(5.2) B = [B tan e L 2 (S,R 2x2 ); B T tan = B tan ). 

As we show below, three main examples of such surfaces are: convex surfaces, 
surfaces of revolution, and developable surfaces without flat regions. 

Lemma 5.6. Assume that S is a simply connected, compact surface of class C 2,1 
with C 1 boundary, and that its shape operator H is strictly positive (or strictly 
negative) definite up to the boundary: 

(5.3) VxeS VteT x S ^|t| 2 < (n(x)r) -r < C\t\ 2 , 

Then S is approximately robust, and more precisely i5.ty) holds. 

Proof. 1. We will prove that every compactly supported, smooth symmetric bilin- 
ear form B t an on S, must be of the form: 

(5.4) Btan = sym Vw, 

for some w € W 1,2 (S, R 3 ). This will clearly imply the Lemma. In |23| . this result 
is proved under an additional assumption that S is closed. The same method, with 
a slight modification, can be applied in our case. For convenience of the reader, we 
present an overview of the argument and for details of calculations we refer to [23] 
and [12], section 9.2. 

Since S is simply connected, it can be parametrized by a single chart r S 
C 2,1 (fi,R 3 ), where f2 C R 2 is a simply connected domain with C 1 boundary. The 
definite form [ffy] j,j:i..2 with gij = dir ■ djr is the pull-back metric on f2 and \/\g\ = 



-y/det[gij] is the associated volume form. Also, the shape operator II expressed in 
the coordinates (x\, X2) G $1 is given by [hij]i t j : i..2> where hij — di(n o r) • djr. The 
inverse of II is denoted II -1 = 2- The mean curvature H of S equals to 

itr ([ 5 y]- 1 n). 

With the above notation, (|5.4|) becomes the following system of partial differen- 
tial equations in fi: 

{dir ■ d\w = Bn 
dir ■ d 2 w + d 2 r ■ d\w = 2Bi 2 
d 2 r ■ d 2 w = B 22 , 

where we set fljj = c^r • B tan djr. Since sym Vw is determined, one concentrates 
on the skew part of Vw. Following [23], we let: 

oj = — -= (diw ■ d 2 r — d 2 w ■ dir) , 
V\9\ 

and we observe that u> must satisfy the equation: 
(5.6) - -^=di i x ~h-u),^) - 2Huj = V{B %3 ). 



151 

The operator V : W 2 ' 2 (fl, R 2x2 ) — > L 2 (Q,R) is a bounded differential operator 
which depends on the geometry of S. The exact expression of D is given in the 
references mentioned before, but for our purposes it is enough to know its stated 
regularity. 
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Now, the following crucial relation between problems (15. 6|) and (|5.5j) is a direct 
consequence of calculations in [23j . 

Proposition 5.7. Assume that \B^\ l& , x ..i € W 2 ' 2 (Q,,R 2x2 ). If \5.6\) has a (weak) 
solution to 6 W 12 (Q, K), f/ien £/ie system i f 5. 5\) has a solution w £ VF 1,2 (f2,R 3 ). 

2. We now show that the hypothesis of Proposition 15.71 is satisfied. Note that 
we have not imposed any boundary conditions on lu, which makes the argument 
easier. Extend first the coefficients and \g\ to and \g\, respectively, defined 
on fl e = {i £ I 2 ; dist < e} for a small e > 0. This extension can be made 

so that [ft 4J ']ij : i..2 satisfies the ellipticity condition (|5.3|) and that h lJ , \g\ and 1/|#| 
stay bounded in Q e . 

In order to prove existence of a solution to (|5.6p . we want to find /o € C^°(O e \0) 
such that the Dirichlet problem 

(5.7) C^=V{B ij ) + f 

has a solution w € W{}' 2 (r2 e , R). The restriction of a> to fl will, clearly, serve our 
purpose. Here the operator C is given: 

is elliptic and self-adjoint with respect to the scalar product: 

(u,C)= f u;CV¥\- 

Therefore, by the classical theory of elliptic operators (see e.g. [7J, section 6.2, 
Theorem 4), (|5 . T[) has a solution if and only if its right hand side satisfies the 
orthogonality condition: 

(P(A J )+/o,C) = 0, 

for all solutions £ € W ' (O e ,M) of the homogeneous problem: = in f2 e . 
The solution space of this problem is finite dimensional, say spanned by a basis 
{Ci) • • • ! Cfc}- For fo 6 C^°(fl e \ fi) consider the functional: 

fc 

L(/ )=^</o,Ci)e i eM fc . 

i=l 

In view of the above, it suffices to prove that L is surjective. 

We now argue by contradiction. Assume that there exists a nonzero a = 
(ati, . . . , at) € M fe orthogonal to the range of L. In other words: 

/ IE a ^i) /ov^I = o V/ ecr(n e \fi), 

which clearly implies that Xh=i a »Ci = in f2 e \ fi. By the Hormander unique- 
ness theorem for second order elliptic equations (see [13], Theorem 2.4), we obtain 
12i=i a id — m ^e, contradicting the linear independence of {(i, . . . , 

In view of Proposition 15.71 this ends the proof. ■ 

Lemma 5.8. Assume that S is rotationally invariant, C 2 up to the boundary, and 
let S have no intersection with its axis of rotation. Then \5. 2\) holds. 
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Proof. 1. After a suitable rigid motion, the surface S can be parametrized by: 

r : (s , si) x [0, 2tt] -» R 3 , r(«, 0) := g(s) 7 (0) + se 3 , 
for a positive function g £ C 2 ([sq, s\], K), e3 = (0,0,1), and 7(0) = (cos 0, sin0, 0). 



As in the proof of Lemma 15.61 we will show that (|5.4p has a solution for £?t a „ 
in an appropriate dense subset of the space in the right hand side of (|5.2p . Given 
w £ W^ 2 (S,R 3 ), write: 

w{s, 0) := o(s, 0)7(0) + 0h'(0) + c(s, 0)e 3 

and also let: 

By = <V • BtandjT. 

The equation (|5.4[) can now be expressed as the following periodic system of 
partial differential equations in (so,Si) x [0, 2tt] (see [27] chapter 12): 

!g'd s a + d s c = B n 
dgb + a = B 22 
g'(d e a-b)+gd s b + dgc = 2B 12 . 

We will prove that (|5.8j) has a solution W 1,2 , periodic in £ [0, 2ir], for Bij being 
finite linear combinations of the Schauder basis for L 2 ([so, s\] x [0, 2tt]) consisting of 
eigenfunctions of Laplacian under the periodic boundary conditions at £ {0, 27r} 
and Neumann boundary conditions in s £ {sq,si}. By density, this will establish 
the lemma. 

2. Differentiating the third equation in s and using the first two equations in 
(15.81) we obtain: 



(5.9) gd 2 b - g"(b + 8 2 b) = 2d s B 12 - d e B 11 - g"d e B 22 =: ^(s, 6). 

Note that ip £ C° and for all s, ij){s, ■) is a finite linear combination of {e }k<N, 
for some integer N independent of s. Hence: 

+00 +00 
&(*, 0) = ^2 b k (s)e tkd and ^(s, 9) = ^ M^Y™ , 

—00 —00 

with ipk = ip-k and tftk = for k > N. Expressing (|5.9p in terms of the Fourier 
coefficients bk and ipk we have: 

b>l- 9 -{l-k*)b k = t±, 
9 9 

Since the coefficients of the above linear equation are continuous in [so,Si], we 
deduce that there exist unique solutions bk £ C 2 ([so, si], E) satisfying bk(so) — 
b' k (so) = 0. Also, bk = 6_fc and bk = for k > N. Concluding, the finite linear 
combination b = J2^k( s )e lke is a W 2 ' 2 solution to (|5.9[) . periodic in 0. One can 
now solve the first two equations in (|5.8|) for a and then for c, obtaining a M^ 1,2 
solution to (15.81) and hence also to (f5T 



Finally, the following result has been proved in [26j . Lemma 3.3: 

Lemma 5.9. Let S be aC 2 developable surface without flat regions. That is, assume 
that for each x £ S the Gauss curvature n(x) = while H(x) 7^ 0. Then S satisfies 
(TOD. 
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6. The recovery sequence: proofs of Theorems 12.21 and 12.3 



In this section, we want to prove Theorems l2.2l and l2.3l that is to define a suitable 
recovery sequence y h . Recall the definition (|2.5[) . With a slight abuse of notation, 
one can write: 

(6.1) Q 2 (x, F tan ) = mia{Q 3 (F ton + c <g> n(x) + rt(x) <g> c); c e R 3 }. 

The unique vector c, for which the above minimum is attained will be called 
c(x,F tan ). By uniqueness, the map c is linear in its second argument. 

Given B tan S B, there exists a sequence of vector fields w h £ W 1,2 (S, R 3 ) such 
that sym Ww h converge in L 2 (S) to B ta n- Without loss of generality, we may 
assume that w h are smooth, and (by possibly reparametrizing the sequence) that: 

(6.2) lim Vh\\w h \\ W 3,~,( S ) = 0. 

Let V e W 2 ' 2 (S,R 3 ) be such that d T V(x) = A(x)t, for all r € T X S, where A € 
14 7l,2 (5, R 3x3 ) is a skew-symmetric matrix field. We approximate V by a sequence 
w h e v^ 2 ^°°(5,]R 3 ) such that, for a sufficiently small, fixed e > 0: 



Yim \\v h -V\\ W 2,2 (S) =0, — — ||w' l |j w/ 2,oo (s) <e , 
(6.3) 71 

lim -/jjieS; « h (a:) ^ War)) = 0. 
h— >o e ft 

The existence of such v' 1 follows by partition of unity and a truncation argument, 
as a special case of the Lusin-type result for Sobolev functions in [TB] (see also 
Proposition 2 in [TP]). 

Define a sequence of rescaled deformations y h <E W 1,2 {S h ° 1 R 3 ): 



y h (x + tn) = x + ^j—v h (x) + Ve^w h (x) 
(6.4) + th/h n(x) + t/h V^(lLv^ an - V(i>' l rf)) (x) 



f2 

th/h V^H T Vw h + th/h V?>d > h {x) + —^hV^d^ h {x). 

/fir. 



Notice that if V £ W 2, °°(S) then one may take v h = V in which case the term 
t/hoV^(Hv^ n - \7(v h n)) is exactly t/h \f^An (see P~3f in Remark g^. 
The vector fields d°' h ,d l ' h £ W hoo (S,R 3 ) are defined so that: 

(6.5) toaVh(\\<P' h \\wi.~w + Wd^Wwi^s)) = 
and: 

lim d 0Jl = 2c (x, B tan - J(A 2 ) tan ) + nA 2 n - \n{n T A 2 n)n in L 2 (S), 

(6.6) \ 2 J 2 

lim d 1Jl = 2c (x, sym (V(Arf) - AII) tem ) + (n T in - n T V(An)) in L 2 (S). 

Lemma 6.1. Assume For £/ie sequence {y h } in \6-4\) the convergences (i), 

(ii) and (Hi) of Theorem \2.2\ hold. 
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Proof, (i) follows by the normalization (|6.2j) . (|6.3j) and ()6.5j) . For (ii) and (iii) notice 
that: 

V h [y h ]=v h + hw h + ^h 2 d 1 ' h 
isym W ,l h/] = -^sym X7v h + sym Vui' 1 + -^-feym Vd 1 ' fi '. 
The proof will be achieved once we establish that: 



(6.7) 



lim y||sym \7v h \\ L 2/ s) = 0. 
h— >o h 



Since the Lipschitz constant of each Vv h is bounded by eo-^rir, anc ^ s y m = 
on the set {x € 5; ^(a;) = V(x)}, we have: 

7p 



|sym Vv h {x)\ < C^=dist (x, {v h = V}) . 



Now we claim that the right hand side above converges to 0, in L°°{S). For oth- 
erwise there would be dist (x h ,{v h — V}) > C^-f^, for some sequence x h g S. 
Consequently, denoting by B x h(r) the ball in R 3 centered at x h and radius r, we 
would obtain: 



fj,{x g S; v h (x) ^ V{x)} > 



SnB x n Qdist (x h , {v h = V}) 



e h 



contradicting ()6.3j) . In the last inequality above we used that the surface S is of 
class C l , with Lipschitz continuous boundary. We thus obtain: 

lim 1 1 sym Vf H^oo^) = 0. 

On the other hand: 

~||sym Wv h \\ L 2 (s) < j^{x g S; v h (x) + V(x)} 1 / 2 ■ ||sym Vv h \\ L - (s) 

The two statements above imply 



Proof of Theorem [2^1 

We will prove that: 

(6.8) lim sup \l h {y h ) < I(V, B tan ) + rj, 

where rj denotes an error quantity, with the property: 

(6.9) 7] -► as e -> 0. 

In view of Theorem l2.1| this will imply (iv) for a recovery sequence obtained through 
a diagonal argument, when eo — * 0. Clearly, the assertions (i) - (iii) will follow as 
well, by Lemma 16. II 
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1. We first look closer at quantities Vh,y h - By Proposition 13. 3\ it follows that: 

(V h y h )(x + tn)n(x) =n+^~ (Man ~ V(v h n)) 

- V7 l n T \7w h + V^d a ^ h + t/hoV^d 1 ^, 
(V h y h )(x + tn)T = Vy h (x + in) (Id + til) (Id + th/hoU^r 

( 6 - 10 ) = |m + Vv h + V^Vw' 1 + th/hoU 

+ t/h VJ l v(llv' t l an - V(v h n)) - th/h ^V{n T Vw h ) 

+ th/h Ve f >S/d°> h + ^hV^Vd 1 ' 11 ) (Id + th/hoUyW, 

for all r G T X S. 

By (|63]1 . ([63]) and one has: \\V h y h ~ Id|| L oo (S h 0) < Ce . It now follows by 
polar decomposition theorem (assuming eo to be sufficiently small), that ^ n y h is 
a product of a proper rotation and the well defined square root of (V 'hy h ) T ^ ' h.y h ■ 
By properties of the energy density function and Taylor expansion, we obtain: 

W{V h y h ) = W f^(V h y h ) T V h y h ^j = W {id + ±K h + 0(\K h \ 2 )j , 

where: 

K h = {V h y h ) T V h y h - Id. 

Clearly: 

(6-11) \\K h \\ L ~ {S H 0) <Ce , 

and so reasoning as in (|3.4p . the above expansion in W yields: 

1 



(6.12) - h WiV h y-) = -Q 3 [^=K- + ^Ofltff) j+^rj- 0(\K^), 
where r\ depends only on ep and satisfies (|6.9p . 



2. Using (|6.10[) we now calculate if 71 . By Error we will cumulatively denote all 
the terms with the property: 



(6.13) lim -^=\\Error\\ L 2 iS h 0) = 0. 
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We start with the tangential minor of K h : 



K£ an (x + tn) = (Id + th/hoU)- 1 



Id + 2— — sym Vv h + 2v / e^sym Vw h + 2th/h U 
h 



+ 2t/h Ve^sym V\Ilv? an - V(v h n)j + -^{Vv h ) T Vv h + t 2 h 2 /h 2 Il 2 



+ 2tVe^/h sym (uVv h ^j + Error 
(Id + th/hoH)- 1 ^ 



{id + th/houy 1 - id 



2sym Vw h + -- r (Vv h ) T Vv h 
br 

+ 2t/h sym V(llvt n -V(v h n) 



2t/h sym (n\7v h ^j 



(Id + th/holl)- 1 + Error, 



where we used the formulae: 



(Id + F) T (Id + F) = Id + 2sym F + F T F, 
Ff 1 FF[~ 1 - Id = F^(F - F^Ff 1 . 



Notice that the quantity Error contains the term ^^-sym ViA, resulting from the 
relaxation of the constraint (|2.2p , (|2.3[) on the small set {v h ^ V}, and other product 
terms, eg: ^-(Vv h ) T V (Jlv^ an — V(v h fi)). The convergence of ^||sym Vv h \\L 2 (s h o) 
to has been proved in (|6.7p . All other terms in Error can be dealt with by 
repeated use of (|6.3p . (|6.2p . Holder and Sobolev inequalities, eg: 



■||(WrV(II^ a „ - V(v h n))\\ L2{s) < C ^-\\Vv h \\ LHs) \\v h \\ w i, Hs) 



< C- 



IL.^II 1 / 2 >,'■ I ' - 



/til 1/2 



<r—\\v h \\ l/2 

SO \\V || W 2,oo( S ) 



as ft -> 0. 



Now, the normal minor of K h is calculated as: 



n T K h (x + tn)n = Ifo t fe a „ - V(v h n) + 2d°- h n + 2t/h d 1 ' h n + Error. 
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The remaining coefficients of the symmetric matrix K (x + tn) are, for r 6 T X S: 



t 1 K h {x + tri)ri = 



h -Vv h + t/h a V^V[Uvt n - V(v h i 



(Id + th/h Tl)- 1 T 



■(ttv? an - V(v h n)Y + {V^d°' h + t/hoV^d 1 

+ ^(nv? an -V(v h n)) T Vv h 
+ t/h V^{uvt n ~V(v h n)) T U 



h\T 



(Id + th/holiy^ + Error 



t/h n T v(llvl n - V(v h n)) + ^(n<„ - V(v h n)) T Vv h 



t/h {uv' t l an ~V(v h n)) T U 



+ {d a > h + t/h a d^ n ) 



l,h\T 



(Id + th/h Iiy X T + Error. 



We leave the estimation in Error to the reader. The convergence of the most 
troublesome term: 

lim i||n T V^ + (IT4 n - V(^n)) T || i2(s , 0) = 0. 

can be proved as in (|6.7p . since the quantity in question vanishes on the set {v h — 
V}. Therefore, \\n T Vv h + (Hv^ an — V (v h n)) T \\ Lao ^ converges to 0, as h — > 0, and 
the displayed convergence follows by the last assertion in (|6.3[) . 

3. In view of (|6.13j) we may now write (with a slight abuse of notation) 

(6.14) ]i m -±—K h = K 1 (x) tan + -^K 2 (x) tan + ((®n + n®<;) in L 2 {S ha ), 

where the symmetric symmetric matrix fields (Ki) ta n € L 2 (S, R 2x2 ) and the vector 
field C G L 2 (S h °,R 3 ) are given by: 

K\{x)tan = B t an — ~ {A 2 ) ta n, 

(6.15) K 2 (x) tan = sym (V(An) - AIl) tan , 

C(x + tn) = c(x, B tan - |(A 2 ) ta „) + — c(x, sym {V(An) - AII) tan ) . 
Further, we observe: 

(6.16) lim -4 I \K h \ A = 0. 

h^o e h J s n 1 1 

Indeed, (16. 14[) implies that -±=K h converges pointwise a.e. in S °. Thus p r |X' l | 4 
converges a.e. to 0. By the boundedness of K h in (|6~TTj) : -^\K h \ 4 < C^\K h \ 2 , 
and the dominated convergence theorem achieves (I6.16p . 

4. Finally, we prove now (|6.8p . By (|6.16p . it follows that the argument of 
Q3 in (|6.12|) converges in L 2 (S h °) to the same limit as 7^K h in ([6~14| . Using 
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Proposition 13. 3\ (I6.14|) and (16. 1| . we obtain: 

r h /2 



i i f r ' 

limsup -rl h {y h ) = limsup-;- / f W{V h y h ) ■ det(Id + th/h Tl) dtdtx 
h^o z h^o e n Js J -ho/2 

< - limsup f Q 3 [ — =K h {x + tft) • det(Id + th/h a U) dtdx 
2 h->0 JsJ-ha/2 \2Ve h ) 

+ Cr) lim sup \ [ \K h \ 2 
h~*o e n J S h 

= ~ I / ^ Q 3 f lim — + Cr? 

2 is i-h /2 2^ / 



< 



-W 

1 r r h °/ 2 / i 

77/7 ^2 a;, ^i (a;) tan + —K 2 (x) tan 

Z JsJ-h /2 V ft 
■h /2 



lim — == 

h->0 2Ve h 



L 2 (S h o) 



dtdx + Crj 



2 Js J -h /2 



Q 2 (x,Ki(x)tan) + 7^-22 (x,K 2 (x) tan ) dtdx + Cr), 



which implies (|6.8p in view of (|6.15|) . 



Remark 6.2. A more careful calculation reveals the exact convergence: 

Yv S n\l h {y h ) = I{V,B tan ), 
h^o e 



for the recovery sequence (|6.4I) . We have used another argument for the sake of a 
more transparent presentation. 

Proof of Theorem [2T3l 

When k — 0, the recovery sequence (for a given V S W /2,2 (S', M 3 ) satisfying (|2.2[) 
and (|2.3[0 is given again by (|6.4p . where we put w' 1 = 0, -Bt an = and n = 0. That 
is: 



y h (a; + in) = x + ——v h (x) + th/h Q n(x) 



+ t/h V^(llvt n ~ V tan {v h n)) (x) + l^hV^d^ix), 



where gJ 1 '" € IT 1,00 ^, R 3 ) satisfies (|6.5p and the second formula in (|6.6p . 

Clearly, y h and V' l [y' 1 '] converge in W 1,2 (S h °) to 7r and V, respectively, as in 
Lemma 16.11 The convergence of the scaled energy follows as in Theorem 12.21 (iv) . 



7. The convergence of minimizers: proof of Theorem 12.51 

Recall that the considered sequence of forces f h € L 2 (S h ,M. 3 ) with zero mean: 
f$h f h = 0, has the form: 

f h (x + tn{x)) = hVe^ det(Id + tll(x))- 1 f(x). 
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Lemma 7.1. Letu h £ W 1,2 (S h ,M. 3 ) be a sequence of deformations such that V h [y h ] 
converges in L 2 (S) to some V : 5 — ► R 3 and let Q h el 3x3 converge to some Q. 
Then: 

} im Ai; f f h -Q h W h -\d) = f f-Qv. 

h^o e" h J S h J s 

Proof. We have: 

4- / f h ■ Q h {u h - id) = i / f H (x) ■ Q h u h {x + tn)~x dtdx 

e 1 h J S h e J s J-h/2 

= \ I f h {x)-Q h ^-V h [y h ] dx= [ f-Q h V h [y h ], 
e Js n Js 

and the result follows. ■ 



Proof of Theorem 12.51 

1. We first show that given any u h £ W 1 ' 2 (S h ,R 3 ) there exists Q h £ 50(3) and 
c h £ R 3 such that w h = (Q h ) T u h - c h - id satisfies: 

(7-1) \\w h f wl ,. 2(sh) < Ch- x I h (y h ). 

Indeed, by Lemma [8~H and properties of the energy density W, it follows that: 



I h {y h ) > Ch- 1 [ dist 2 (Vu h , 50(3)) > Ch- 1 f \Vu h - R h Tr\ 2 
Js h Js h 

(7.2) >Ch- x \ \(Q h ) T Vu h -Id| 2 - Ch' 1 f \{Q h ) T R h Tr-ld\ 2 

Js h Js h 

> Ch- 1 f \Vw h \ 2 - Ch- 2 I h {y h ). 
Js h 

Actually, the assumption of smallness of h~ 3 E(u h , S h ) cannot be expected to hold 
here. In this general case one exchanges the 50(3)-valued matrix field R h with 
-R' 1 £ W 1 ' 2 (S, M 3x3 ) given in the proof of Lemma 15711 Then Q h for which the above 
estimates are true may be taken as a rotation in 50(3) with minimal distance from 

By (f72| it follows that: 

||V^||| 2Csh) < Chl h {y h ) + Ch- 1 I h (y h ), 

which implies (|7.ip in view of the Poincare inequality, for an appropriately chosen 
constant c . A proof of the uniform Poincare inequality on S h can be found, for 
example, in [PT] . 



2. Notice that by the definition of m h we have: 

rh/2 
>S J-h/2 



tI f h {z)-Q h z dz = h\fe T ( / 1 f(x) ■Q h (x + tn) 

< l Js h J S J-h/2 



< m h . 
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Therefore, in view of (12.71) and (17.11) we obtain: 



C h u h 



J h (y h )-I h {y h ) = m h - T / / 



h 



= ~f (Q h ff h ■w h + m h -U f h - Q h z dz 

> -r / (Q h ff h ■ ™ h > -C/i 1/2 v^||/|U a(s) ||^|| L2(sh) 
n Js h 



(7.3) 



> -CV^I h (y h ) 1/2 . 

We now prove the first claim of the theorem. Taking u h (z) — Qz for any Q £ A4, 
we notice that J h (y h ) — 0. Hence inf J h < 0. The lower bound of \J h follows 
from ([73]) : 

1 /2 

(7.4) ^)>^)-c(^)) , 
which proves (i). 

3. To prove (ii), let u h be a minimizing sequence of jn-J h , as defined in (|2.9[) . 
Then J h (y h )} is bounded, and therefore, by (|7.4j) {pr-/ 7l (y' 1 )} is also bounded. 
The convergences of y h , V h [y h ] and ^sym VV h [y h ] follow from Theorem O In 
particular: 

(7.5) Um inf \l h {y h ) > I(V, B tan ). 

h-^Q e' 1 

The strong convergence of ^sym W h [y h ] is deduced from the strong convergence 
of the sequence symG^ in Lemma 13.61 This last result is in turn implied by 
the convergence of J s Qy,{G h ) (valid because the sequence is minimizing), positive 
dcfiniteness of Q3 on symmetric matrices, and the weak convergence of G h . Since 
the details are exactly the same as in [10] section 7.2, we omit them. 

We now prove that the limit Q of any converging subsequence of Q h belongs to 
M. By (|2"77|) we have: 

(7.6) = (m- f -f 1 f{x) ■ Q h u h dtdx] 

JS J-h/2 I 



f f h ■ Q h (u h -id) + -^=( m - f f ■ Q h x da J 
Js h ve h V Js J 



The first term above is bounded, as it in fact converges to — J s f-QV, by Lemma l7TTl 
The quantity is brackets in the second term converges to m — J s f ■ Qx. Therefore, 
if Q £ M, this last quantity is uniformly positive, and the second term above 
converges to +00 (as h/Ve^ — > 00). We observe that, in this situation, ^hJ h (y h ) 
must converge to +00, contradicting (i) and thus proving that Q e M. 
In view of (|7.5[) . (|7.6[) also implies: 

hminf \j h {y h ) > hminf \l h {y h ) - f f • QV > J(V, B tan ,Q). 
h^a e h— >o e J s 
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The fact that the limit (V, Btan,Q) minimizes the functional J is now a standard 
consequence of the above inequality. Indeed, if: 

J(V,B tan ,Q) < J(V,B tan ,Q)-e 

for some V £ W 2 - 2 (S, R 3 ) satisfying some B tan £ B, Q £ M and e > 0, 

then for a related recovery sequence y h there would be: 

lim \j h {Qy h ) = J(V,Btan,Q) < J{V,B tan ,Q) - e < liminf \j h {y h ) - e, 

which contradicts (|2.9p . 

Finally, (iii) follows exactly as (i) and (ii). ■ 

Remark 7.2. 1. A dead load (versus a "live load") is any external force which only 
depends on the reference configuration point, and not on the deformation itself. An 
important feature of dead loads, discussed first in [21], is the following. If the load 
is in a certain average sense compressive, it is advantageous for the body to perform 
a large rotation rather than undergo a compression. Our analysis identifies A4 as 
the set of candidates for such rotations, which are expected to minimize the total 
energy J h among all rigid motions of the body. 

This phenomenon may happen even if the average torque of the force is zero: 

(7.7) ( f{x) xidi = 0. 

Js 

Note that vanishing of the average torque is necessary for Id £ M. , since (|T. T[) can 
be written as: J s f ■ Fx = for all F £ so(3). However, it is not sufficient, and 
if Id ^ M. then we observe that the minimizers of J h will not be close to Id. In 
general, the body chooses an infinitesimal isometric displacement V and a rotation 
Q £ M. which is energetically advantageous in response to the force /. That is, 
those rotations which allow for a better alignment of infinitesimal isometries with 
the direction of the dead load, are preferred. 

2. The assumption on the sequence of forces f h can of course be weakened. For 
example, consider f h (x + tn) = dct(Id + tH(x))^ 1 f h (x) and let h ^rr: f h converge 
weakly in L 2 (S) to some / £ L 2 (S, M 3 ). In this situation, one needs to enforce extra 
assumptions on the asymptotic behavior of the maximizers of the linear functions 
50(3) 3 Q i— > J s f h ■ Qx dx with respect to M, to exclude certain degenerate 
cases. The analysis is as in the proof of Theorem 12.51 and we leave the details to a 
courageous reader. 

3. The lower bound on J and existence of its minimizers can be proved inde- 
pendently, and under the following weaker assumptions: 

[f = and / f(x) ■ QFx dx = VQ £ M V-F £ so(3), 
Js Js 

which can be seen as the linearization of (|2.8p . although it makes perfect sense for 
any closed nonempty subset M. C 5*0(3). Indeed, the second equality above follows 
by differentiating the expression J s f(x) ■ Qx dx at Q £ M. and using that so(3) 
is the tangent space to 50(3) at Id. We present the proof of coercivity and the 
attainment of the minimum by J and J under this condition, for arbitrary M, in 
Appendix C. 
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8. Appendix A - an approximation theorem on surfaces 

For a given vector mapping u £ W x,2 (il, M. n ) denned on an open subset ft C R n , 
denote: 

E(u,to)= [ dist 2 (Vu(x),S'0(3)) dx. 
Jn 

Lemma 8.1. Let u £ W 1 ' 2 (S h ,R n ) and assume that h~ 3 E(u, S h ) is sufficiently 
small. Then there exists a matrix field R £ W ' (S, R 3x3 ), such that: 

R(x) e SO(3) Vx e S, 

and a matrix Q £ 5*0(3) with the following properties: 

(i) \\Vu-Rn\\ L2{sh) <CE(u,S h y/ 2 , 

(ii) \\VR\\ LHS ) <Ch-^ 2 E(u,S h f/ 2 , 

(iii) \\Q T R - U\\ LP{S ) < Ch- 3 / 2 E(u, S h )^ 2 , for all p £ [1, oo), 
where C is independent of u and h (but may depend on p). 



The proof of Lemma 18.11 uses the following nonlinear quantitative rigidity esti- 
mate by Friesecke, James and Miiller: 

Theorem 8.2. |9J Let VL C M™ be an open, bounded domain with Lipschitz bound- 
ary. Then, for every u € W ' (fi,R n ) one has: 



[ \Vu(x) - R\ 2 dx < CE(u,Cl), 
Jn 



mm 

Reso(n) J Q 

where the constant C depends only on fl. In particular, C is invariant under 
dilations of Q, and it is also uniform for the uniform bilipschitz images of a unit 
ball m R". 

Proof of Lemma 18.11 

1. For x £ S define 'balls' in S and the corresponding 'cylinders' in S . 

D x , h = B{x, h) n S, B x , h = it-^D^h) n S h . 

The main observation is that Theorem l8.2l mav be applied on each set B x .h, yielding 
matrices R x ,h € 50(3) such that: 



(8.1) 



/ \Vu{z)-R x j l \ 2 dz<CE{u,B x ^ h ), 



with uniform constant C (independent of x or h). 

2. Let $ £ C£°([0, 1)) be a nonnegative cut-off function, equal to a nonzero 
constant in a neighborhood of 0. For each x £ S define the function n x : S h — > R: 

6{\irz-x\/h) 
Vx[Z> f sh 4(\*y-x\/h)dy- 
Then rj x (z) = for z B x ^ and: 



(8.2) / r) x (z)dz = l, <Ch- 3 , \\V x r) x \\ L °° <Ch~ 



'/•' 1 ■ ' ' 1 - — '• • !l 7.' II L ' •• ■' ' 

IS h 

The last inequalities follow from the lipschitzianity of dS. In particular, the de- 
nominator function in the definition of r\ x has Lipschitz constant Ch 2 , and hence: 

V, [ 0(\ny - x\/h) dyj < Ch-\ 
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Consider the matrix field R £ W 1 ' 2 (S,R Sxi ) 



R(x) 



r] x {z)Vu{z) dz. 



By the first two statements in (J872J) we obtain: 



(8.3) \R(x) - R x , h \ 2 = 
Similarly: 



rj x (z)(\7u(z) - R x . h ) dz 



<Ch- A E(u,B xA ). 



|Vi?0r)| 2 = 



(8.4) 



('V x r) x )'Vu 



S 1 ' 



(V x r] x )(yu - R Xth ) 



<[ I |V«- R x ,hf <Ch- 5 E(u : B xJl ) : 

and likewise, for any x l £ D x ^: 

(8.5) |V7?(x')| 2 < Ch- 5 E(u, 2B x<h ) 

with 2B xJl = Tr ^(D X 2h) n S h . Therefore, in view of the lipschitzianity of OS and 
by the fundamental theorem of calculus: 

(8.6) \R(x") - R(x)\ 2 < Ch- 3 E(u, 2B x<h ) Vx" e D x<h . 
Combining 4EH]) with (JO]) and (8j| yields: 

\Vu(z) -Rtt{z)\ 2 dz 



(8.7) 



< 2 



\Vu-R xM \ 2 + 



\R(x)-R x , h \ 



\Rir - R(x) 



< CE(u,2B x>h ). 



Now cover 5 by {-Dzj.fcJiJi so ^ na ^ the coverm g number of the family {2B Xu h}^i 
is independent of h. Summing the inequalities in (|8.7[) over i = 1 . . . N proves: 

\Vu-Rtt\ 2 < CE(u,S h ). 

Also, integrating (|8.5p over D Xit h and summing over i = 1 ... N gives: 

(8.9) f \S7R\ 2 < Ch- 3 E(u 7 S h ). 

Js 

3. Notice that by (fOj) . for every x G 5: 

dist 2 (R(x),SO(3)) < Ch- 3 E(u,S h ). 
Hence, if E(u, S h )/h 3 is sufficiently small, we may define: 

R(x)=V so(3) (R(x)) 

where Vso(3) is the orthogonal projection onto the compact manifold 50(3). Clearly 
R : 5 — > 50(3) is a W 1 ' 2 matrix field and since: 

\R(x) - R(x)\ = dist (£(x), 50(3)), 
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by (18. 8|) we conclude that: 

/ \Vu~ Rir\ 2 <C ( [ \Vu-Rir\ 2 +[ dist 2 (Vu, 50(3)) ) < CE(u, S h ), 
Js h \Js h Js h J 

which proves (i) in Lemma |8. 11 The bound (ii) is deduced directly from 



4. To deduce (iii), define first the intermediate matrix Q as the average of R on 
S. Using the Sobolev and Poincare inequalities, together with (ii) we obtain, for 
every p > 2: 

(8.10) (J \R - Q\A ^ <C\\R- Q\\ 2 W ^ [S) <C J | Vi?| 2 < Ch~ 3 E(u, S h ). 

Now, take Q S SO(3) such that \Q - Q\ = dist (Q,SO(3)). As before, (18T01) 
remains true if Q is replaced with Q. Clearly, the same bound must also hold for 
p G [1,2), and so we conclude that: 

Vpe [l,oo) \\R-Q\\Ip {s) < Ch- 3 E{u,S h ). 

The above easily implies (iii). H 



9. Appendix B - the T-convergence setting 

We first recall the notion of T-convergence of a sequence of functionals T h : 
X — ► M, defined on a metric space X. Namely, T h T-converge, as h — * 0, to some 
T : X — ► K provided that the following two conditions hold: 

(i) For any converging sequence {x h } in X one has: 

(9.1) T\ lima;' 1 ) < liminf T h (x h ). 

\h— >0 J h— >0 

(ii) For every x £ X , there exists a sequence {x h } converging to x, such that: 

(9.2) T{x) = \imT h {x h ). 

When X is only a topological space, the definition of T-convergence involves, nat- 
urally, systems of neighborhoods rather than sequences. However, when the func- 
tionals T h are equi-coercive and X is a reflexive Banach space equipped with weak 
topology, one can still use (i) and (ii) above (for weakly converging sequences), as 
an equivalent version of this definition. For details, we refer the reader to [B]. 

Proof of Corollary l2~4l 

We only prove (i) , in the case when the product space in the domain of T is equipped 
with the strong topology. The other statements follow the same. 

To obtain (|9.1[) . we take a sequence oiW 12 (S ho ) vector mappings {y h } such that, 
writing B^ an = isym VV h [y h ], the sequence {T h {y h , V h [y h ], B? an )} is bounded, 
and such that y h , V h [y h ] and B^ an converge to some y, V and B tan (in W 1 > 2 (S h °), 
W 1 ' 2 (S) and L 2 (S) respectively). By Theorem 12.11 we obtain a sequence of nor- 
malized deformations y h — (Q h ) T y h — c , converging to w. Moreover, U /l [y' 1 ] and 
^sym W h [y h ] converge to V and weakly to B tan , respectively. Notice now that: 

\Q h -Id| = Ch- 1 Ve^\\VV h [y h ] ~ {Q h ) T VV h [y h } \\ L2 < Ch^Ve*. 
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In particular, Lemma 13. II remains true if we put Q h = Id, for all h. Consequently, 
all the assertions of Theorem l2.1l still hold for y h = y h — c h (possibly after modifying 
the constants c ). 

Now, V A,l [?/' 1 ] — V h [y h ] — h/Ve^c h is bounded, so c h converge to 0, as h — > 0. 
On the other hand c h = y h — y h converge to y — n. Hence y = tt. Moreover 
VV h [y h ] = VV h [y h ], so W = VV and B tan = B tan . By Theorem O (iv) we 
conclude that: 

Hy, V, B tan ) < limmf T\y h , V h , B* n ) 

h— vu 

which proves (|9.1|) . 

The second requirement for T-convergence (19.21) follows directly from Theorem 
O in view of (Ol. ■ 



Wc remark that in presence of external forces, the results of Theorem 12.51 can 
also be formulated in the language of T-convergence, similarly as above. 

10. Appendix C - on coercivity of the generalized von Karman 

FUNCTION ALS J AND J 

In this section, we consider the functionals: 

J(V, B tan ,Q) =~ J Q 2 (x, B tan - |(A 2 ) ta „) 

+ li Is 22 {x ' (V(j4 " } ~ AU)tan) ~ l s f ' QVl 
J(V, Q) ~ j s Q 2 (x, (y(Afi) - AU) tan ) - J f.QV, 

defined for infinitesimal isometries V, matrix fields B ta n £ B and rotations Q G jM, 
where A"! is an arbitrary closed and nonempty subset of SO (3). We prove that J and 
J attain their finite lower bounds under the following assumptions on / G L 2 (S, R 3 ): 

(10.1) // = and [ f{x) ■ QFx dx = VQ G M VF G so(3). 

As mentioned in Remark 17.21 the second condition above is a consequence and 
linearization of (|2.8|) , with A4 defined by that formula. 



Lemma 10.1. Assume that S is of class C 2,1 . Then for every V G W 2 ' 2 (S, R 3 ) 
satisfying 12. 2j) and S2.3\) there exist D G so(3) anrf d G R 3 , so that: 

\\V - (Dx + d)\\ 2 w2 , 2{S) <C j \(V(An) - AU) tan \ 2 . 

J s 

Proof. 1. We first prove that J s |(V(An) — AII) ta „| 2 = implies for a W 2 ' 2 in- 
finitesimal isometry V to have the form V(x) — Dx + d, with D G so(3) and 
de R 3 . 

To see this, let c G W 1 ' 2 {S,R 3 ) be such that: 

A(x)t = c(x) xr Vx G S Vt g T x 5. 

Since A represents a gradient, it follows that d T c x 77 = c^c x r for all r, 77 G 1^5. 
In particular, for any r and 77 such that r x 77 = n, one has: 

(<9 T c) • n = —(d T c x 77) -t = ~(d v c x r) ■ r = 0. 



34 MARTA LEWICKA, MARIA GIOVANNA MORA AND MOHAMMAD REZA PAKZAD 



On the other hand: 

= (d T {An) - AUr) tan = (d T (c x n) - AUr) tan = (d T c) x ft. 

Hence d T c = on S, which yields the claim. 

2. We prove the result. Arguing by contradiction, we assume that for a sequence 
of infinitesimal isometries V h G W 2 ' 2 (S, R 3 ) there holds: 

distvy2,2 (s) (v h , {Dx + d- De so(3), d G R 3 }) = 1 

(10 - 2) r, „ ,2 

and lim / \{V A h )n\ = 0. 

Since the second condition above involves only higher derivatives of V , we may 
without loss of generality replace the first condition by: 

(10.3) 11^11^2,2(5) = 1 and {V h ,Dx + d) w2 2{s) =0 WD G so(3), d G M 3 . 

In particular, converges weakly in iy 2 ' 2 (S') (up to a subsequence, which we do 
not relabel) to some vector field V, still satisfying (|2.2p and (|2.3p . By (110.2)) and 
the weak lower semicontinuity of the L norm, we deduce that J s \(VA)n\ = 0. 
Hence, in view of the first part of the proof and the second condition in (|10.3[) . it 
follows that V = and so: 

(10.4) lim II^H^s) =0. 
By the estimate (|4.4I) and (|10.4p we may deduce: 

(10-5) lim||^JIw2,2 (S ) =0, 

where V^ n = V h — (V h n)n. Observe that: 

/ \(VA h )n\ 2 = [ V {YWtn - V{V h n)) - A h Il " 
Js Js 

= [ V 2 (V h n) + (A h U - UA h ) tan - (VU)V t h an + (V h n)U 
Js 

Therefore: 

\\V 2 (V h n)\\ LHS) < c(\\(VA h )n\\ L2(s) + \\V h \\ w ^ {s} ) , 

and in view of (|10.4|) and the assumption (|10.2[) we also get: 

lim \\V h n\\ W 2,2 (s) = 0. 

The above together with (|10.5|) contradicts (|10.3|) and proves the lemma. H 

Lemma 10.2. Assume H10.1\) and let S be of class C 21 . Then the junctionals J 
and J, defined for V £ W 2 > 2 {S,R 3 ) satisfying fSJfy . fOj) . and B tan eB, Q G M, 
are bounded from below and attain their infima. 

Proof. 1. Let V G W 2 ' 2 (S, R 3 ) be an infinitesimal isometry. By Lemma Il0.1[ 
positive definiteness of Q2 (on symmetric matrices) and (|10.1j) . we obtain: 



2 



nnc , J(V) > C||F||^ 2 , 2(S) - f.QV= C\\V\\ 2 W2 , 2{S) - f-QV 

(10.6) Js 

>C\\nw^s)-\\f\\LHS)-\\V\\LHS), 
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for an appropriate modification V = V — (Dx + d) . Hence the lower bound on J 
(and J) follows. 

2. Let now (V h , B^ an , Q h ) be a minimizing sequence of J. Clearly, a subsequence 
of Q h converges to some Q E M.. 

Using (|10.6p and applying the positive definiteness of Q2 to the first term in J, 
there follows the (uniform in h) boundedness of the following expressions: 

(10.7) (c\\V h fw*.HS) - \\fh*(S) ' \\V h \\ms } ) + C 

Again, we put V h = V — (D h x + d h ) and apply Lemma TIP, f I In particular, the 
sequence V h is bounded in W 2 ' 2 (S) and so it converges (up to a subsequence), 
weakly in W 2,2 (S), to an infinitesimal isometry V. Further, the matrix fields A h = 
A h - D h converge weakly in W 1 ^ 2 {S) to the field A satisfying and ([23| . 
Notice that: 

(A h ) 2 = (A h ) 2 + {D h f + (D h A h + A h D h ). 

Hence the boundedness of the second term in (|10.7|) results in the L 2 (S) bounded- 
ness of: 

BL ~ \ ((D h ) 2 + (D h A h + A h D h )) tan = Bt n - |sym V {{D h ) 2 x + 2D h V h (x)) . 

We may now conclude that a subsequence of the above sequence of symmetric ma- 
trix fields converges, weakly in L 2 (S), to some B tan £ Thus, B' t l an — ^((A h ) 2 ) tan 
converges to B tan - § (A 2 ) tan . 

By the weak lower semicontinuity of both quadratic terms in J we conclude that 
J(V, BtanQ) realizes the infimum of J. Likewise, J(V, Q) realizes the infimum of 
J, had V h been a minimizing sequence of J. M 
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